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Abstract. Dark matter haloes from cosmological N-body simulations typically have triaxial 
shapes and anisotropic velocity distributions. Recently it has been shown that the velocity 
anisotropy, /3, of cosmological haloes and major merger remnants depends on direction in 
such a way that /3 is largest along the major axis and smallest along the minor axis. In 
this work we use a wide range of non-cosmological N-body simulations to examine halo 
shapes and direction-dependence of velocity anisotropy profiles. For each of our simulated 
haloes we define 48 cones pointing in different directions, and from the particles inside each 
cone we compute velocity anisotropy profiles. We find that elongated haloes can have very 
distinct velocity anisotropics. We group the behaviour of haloes into three different categories, 
that range from spherically symmetric profiles to a much more complex behaviour, where 
significant differences are found for (3 along the major and minor axes. We encourage future 
studies of velocity anisotropics in haloes from cosmological simulations to calculate /3-profiles 
in cones, since it reveals information, which is hidden from a spherically averaged profile. 
Finally, we show that spherically averaged profiles often obey a linear relation between (3 and 
the logarithmic density slope in the inner parts of haloes, but this relation is not necessarily 
obeyed, when properties are calculated in cones. 

Keywords: galaxy dynamics, dark matter simulations, dark matter theory 



Contents 



1 Introduction 2 

2 Simulations 3 

2.1 Simulation code 3 

2.2 Simulation I - An instantaneous change in the potential 4 

2.3 Simulation II - Exchange of energy 5 

2.4 Simulation III - A cold collapse experiment with substructure 5 

2.5 Simulation IV - Major mergers 6 

2.6 Simulation V - Unstable Osipkov-Merritt models 6 

2.7 Simulation VI - A halo and its subhaloes 6 

3 Analysis methods 7 

3.1 Analysing particles in cones 7 

3.2 Determination of halo shapes 7 

4 The behaviour of velocity anisotropy and halo shape 7 

4.1 Category 1: Haloes with well-behaved /^-profiles 7 

4.2 Category 2: Haloes with complex /3-profiles 8 

4.3 Category 3: The spherically symmetric haloes 8 

4.4 A summary of the categories 9 

5 Relating the results to cosmological haloes 10 

5.1 Comparison with a cosmological halo 10 

5.2 Spatial anisotropy of galaxy kinematics in galaxy clusters 11 

6 Spherically averaged profiles 11 

6.1 Velocity anisotropy profiles 11 

6.2 7-profiles 11 

6.3 Halo shapes 12 

7 Comparison with the attractor 15 

7.1 Motivation for the attractor 15 

7.2 Comparing haloes with the attractor 15 

7.2.1 Different projections 16 

7.2.2 Comparison with cosmological pseudo-phase-space densities 16 

7.2.3 Energy exchange in the relaxation processes 16 

7.3 An overview: which structures are on the attractor? 19 

7.4 The relevance of the attractor 19 

8 Summary 19 



- 1 - 



1 Introduction 



The gravitational potential of galaxies and galaxy clusters are dominated by dark matter 
haloes. Such haloes have only been observed through their gravitational impact on other 
objects: e.g. in the Bullet Cluster [1, 2], the motion of galaxies in clusters [3] and in the 
cosmic microwave background [4]. Therefore, the most fruitful way to study the detailed 
dynamics of dark matter structures is through numerical simulations. 

The formation and evolution of the large-scale structure in the universe can be simu- 
lated with N-body codes, where the density field and velocity field of the dark matter are 
represented by coUisionless particles. In such simulations it is found, that the dark matter 
haloes follow a universal density profile, usually parametrised as an Einasto-profile [5] or a 
NFW-profile [6]. A different universality is the pseudo-phase-space density profile, p/crf^^, 
which is a power law in radius {p is the density and o"rad is the radial velocity dispersion) 
[7-10]. It is worth noting that particles in subhaloes and streams do not follow this power 
law behaviour [11]. Subhaloes are also responsible for large peaks in the true 6-dimensional 
phase-space density distributions in haloes from cosmological simulations [12-16]. 

The mass distributions of haloes from cosmological simulations are triaxial ellipsoids 
with a major axis (with length a), an intermediate axis (with length b) and a minor axis 
(with length c). If a ~ 6 > c a halo is oblate {pancake shaped), if a > 6 ~ c a halo is prolate 
{cigar shaped), and if a ~ 6 ~ c a halo is spherical. High-redshift haloes are typically more 
elongated than those at low redshift, and it is also a trend that light haloes are closer to 
being spherical than massive haloes [17]. One possible mechanism that likely is important 
for creating elongated haloes is mergers, that can produce both prolate and oblate haloes 



Dark matter systems can have anisotropic velocity distributions. The anisotropy is 
characterised by the velocity anisotropy parameter, 



where (Trad is the velocity dispersion in the radial direction, and Utan is the total tangential 
velocity dispersion. /3 is positive in regions with radially dominated anisotropy, and negative 
in regions with tangentially dominated anisotropy. The inner parts of haloes in cosmological 
simulations have /3-profiles (calculated in spherically averaged bins) increasing from in the 
deep interior to about 0.3 at the core radius [11, 20-22]. The /3-profiles in the outer parts 
exhibit different behaviour from halo to halo. 

From the outcome of a wide range of simulations [23], a relation has been suggested 
between the spherically averaged value of /3 at a given radius (r) and the slope of the density 
profile. 



The relation, /3(7) = —0.2 x (7 + 0.8), is valid in the inner parts of haloes in cosmological 
simulations [24]. A more recent relation, the attractor, was found in several different struc- 
tures, which were exposed to an artificial process that instantaneously interchanged energy 
between particles [25] (HJS). The role of the attractor in cosmological simulations, however, 
remains undetermined. 



[18, 19]. 
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Table 1. The softening, the total mass and the gravitational constant 
all the initial structures. 



in the simulations. Vg = 1 in 





Description 


Softening 


Mass 


I 


G- variations 


0.0050 


1 


II 


Energy exchange (HJS) 


0.0050 


1 


III 


Cold collapse 


0.0025 


1 


IV 


Multiple mergers 


0.0050 


128 


V 


Unstable models 


0.0050 


1 


VI 


Substructure 


0.0015 


1 



There exists evidence that /3-profiles are not spherically symmetric functions (or func- 
tions that are constant along density or potential contours, if the structure is non-spherical) 
in all dark matter haloes. In a detailed study of remnants of collisionless mergers, it has 
been shown that the velocity anisotropy profiles of merger remnants exhibit different be- 
haviour depending on whether they are calculated along the major or the minor axis [26]. 
The behaviour of /3 were clearly correlated with the shape and the orientation of the haloes, 
in such a way that the largest velocity anisotropy was found along the major axis. Similar 
asymmetric /3-profiles are also found in cosmological simulations [27]. 

The aim of this article is to study the behaviour of direction-dependent velocity anisotropy 
profiles and halo shapes in a set of non-cosmological simulations. In section 2 we describe our 
simulations, and in section 3 we describe our analysis methods. Next (section 4) we study the 
direction-dependence of /5-profiles together with halo shapes, and we compare our results to 
cosmological haloes in section 5. In section 6 we report halo shapes and spherically averaged 
/3- and 7-profiles. Finally we compare haloes with the attractor from HJS (section 7). 

2 Simulations 

The overall purpose of our simulations is to expose haloes to a range of perturbations, so we 
can study velocity anisotropy profiles and halo shapes in haloes formed in several different 
ways. Simulation I is an artificial simulation, that involves an instantaneous change in the 
gravitational potential, simulation II is the energy exchange perturbation from HJS, simula- 
tion III is a collapse experiment, simulation IV involves multiple major mergers, simulation 
V involves unstable haloes and simulation VI shows the effect of substructure in a halo. An 
overview of parameters is given in Table 1, and an overview of the initial structures is given 
in Table 2. 

2.1 Simulation code 

To run our simulations we used the N-body simulation code, Gadget-2 [28, 29]. The simula- 
tions were run in a non-cosmological Newtonian box, and we only used collisionless particles. 
For all the simulations the spline softening, r], implemented in Gadget-2 was used. The 
time-step of each particle was calculated as At = (2r/e/|a|)^/^ (for a discussion of time-step 
criteria see [30]), where |a| is the magnitude of the acceleration and e is the accuracy param- 
eter, which was set to 0.05 in our simulations. In all simulations the energy conservation was 
better than 1.0 %. 



-3- 



Table 2. An overview of the structures in the simulations. The Colour and Symbol columns refer to 
how structures are presented in figure 7-13. 



Same collision axis in the two last simulations. 



Name 


p 


P 


Sim. 


Colour 




Symbol 


A 


Hernquist 





I and II 


maroon 







B 


1/(1 + rf 





I and II 


red 




o 


c 


Hernquist 


r^/ri + r^), ran = 1-0 


I and II 


orange 




< 


D 


1/(1 + r)3-5 


r^/ri.S^ +r2), r.n = 1.5 


I and II 


blue 




X 


E 


Hernquist 


0, Gaussian 


I and II 


pink 




* 


F 


Hernquist 


1/2, Gaussian 


I and II 


green 




+ 


G 


Hernquist 


r^/(l + r^), Gaussian 


I and II 


black 







H 


Hernquist 


1 — r'^/ (1 + ^^^l, Gaussian 


I and II 


yellow 




V 


I 


Hernquist 




HI 


black 




o 


J 


1/(1 + r)5 


- 


HI 


black 




X 


K 


Hernquist 





IV 


red 







L 


1/(1 + r)5 





IV 


red 




* 


M 


1/(1 + r)3 





IV 


red 




< 


N 


Hernquist 





IV^ 


blue 




o 





1/(1 +r)5 





IV^ 


blue 




* 


P 


1/(1 +r)3 





IV^ 


blue 




< 


Q 


Hernquist 


rV(0.2^ + r^), ran = 0.2 


V 


green, IC: 


blue 


<l 


R 


1/(1 +r)5 


rV(0.22+r2), ran = 0.2 


V 


green, IC: 


blue 


* 


S 


Hernquist 





VI 


red, IC: black 


0, IC: o 



2.2 Simulation I — An instantaneous change in the potential 

The first of our simulations will show how structures respond to a process, where the gravita- 
tional potential changes instantaneously. Such a change perturbs the accelerations in contrast 
to the perturbations of the velocities performed in HJS (see section 2.3 for details). 

First eight structures with different velocity anisotropics and density profiles were cre- 
ated. The density profiles were of the form, 

Pir) = (2.1) 

(r/r,)-«(l+r/r,)-« 

where rg is the scale radius and po is a normalization constant determining the total mass. 
Six structures followed a Hernquist density profile (C = 1 and C = 3) [31], one structure had 
(^, C) = (0,3.5) and one structure had (^,C) = (0,5). 

Two of the structures were created using Eddington's formula for /3 = [32], and two 
structures were Osipkov-Merritt models [33, 34] with velocity anisotropy profiles given by, 

/5W = (2.2) 

' an ~r ' 

where ran is the anisotropy radius. 

For four of the Hernquist structures the initial velocity distributions were Gaussian 
distributions with velocity dispersions, a^nd and <Ttan, calculated from the Jeans equation [35]. 
Structures with initial = 0, Anitial = 1/2, Anitial = '^^/(l + r^) and initial = 1 - r2/(l + r^) 
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were created. Structures generated with this method are not in perfect equihbrium, so 
simulations were run for 100 time units, where the structures had time to equihbrate. We 
also ran such test simulations for the structures generated with Eddington's method and for 
the Osipkov-Merritt models to assure equilibration. 

All structures had rg = 1, and po was chosen such that the total mass inside r = ISrg 
was 1. When structures are created in this way the dynamical time, 1/y/Gp, for particles 
inside r = 13rs is smaller than 100 time units, which is the duration of all our simulations. 
The initial conditions were generated with G = 1. The structures used in this simulation are 
identical to those used in HJS. The initial density profiles and velocity anisotropy profiles 
are summarised in table 2. 

The next step in our setup was to run a simulation for 100 time units, where the 
gravitational constant was set to G = 0.8. In this time the structures expanded due to the 
lowering of the potential. Next the gravitational potential was increased by setting G = 1.2, 
which caused a contraction in the following 100 time units. We kept changing between 
G = 1.2 and G = 0.8 every 100 time units until 2000 time units had passed. In the last of 
these simulation we had G = 1.2. An additional simulation with G = 1.2 was run for 100 
time units, so the structures had more time equilibrate. 

2.3 Simulation II — Exchange of energy 

Simulation II is the simulation from HJS, where an exchange of kinetic energy between 
particles in the same spherical bin was performed. The experiment was performed for the 
haloes described in section 2.2. In each perturbation the kinetic energy of each particle was 
multiplied by a uniformly chosen random number in the range [0.25, 1.75] . Conservation of the 
total energy was taken care of by scaling all the kinetic energies by a constant (the constant 
was typically between 0.98 and 1.02). After this perturbation the haloes were evolved in time 
with a N-body simulation code. 20 of such perturbations (followed by time evolution of the 
system) were performed. See more details in HJS. 

2.4 Simulation III — A cold collapse experiment with substructure 

In this simulation we want to assemble a structure through a violent relaxation process 
(similar to [36]), that mimics the way structures are assembled when they collapse in the 
early universe. Our setup will be a main halo, which contains several compact substructures. 
Initially, all particles had a velocity of zero, so this simulations is therefore effectively a 
collapse simulation, where the substructures break the spherical symmetry. 

We distributed 5 x 10^ particles according to a Hernquist density profile with = 1 
and a cutoff at 200rs. Next 5 x 10^ particles, with the same total mass as the main halo 
particles, were distributed in 24 identical subhaloes, that also followed Hernquist profiles, 
but with a scale radius of 0.5 and a cutoff radius of 5. The positions of the subhaloes were 
sampled in the same way as the particles in the main halo. The total mass in the simulation 
was 1. With this setup the time scale of the collapse is 2.45 times larger for the main halo 
than in the subhaloes. 

We ran the simulation for 200 time units, which corresponds to 200 dynamical times at 
the scale radius for the initial structure. An additional simulation was run, where the initial 
density profiles of the main halo and the subhaloes were. 



p{r) 



Po 



(2.3) 
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We defined rg = 1, and po was found by defining the total mass of the halo to be 1. 

Note that this simulation is a physically realistic experiment, in the sense that no 
artificial perturbations has entered the simulation, even though the initial condition for the 
collapse is different from that of cosmological simulations, e.g. because the initial conditions 
have a cuspy density profile. 

2.5 Simulation IV — Major mergers 

Now we will build up a structure through major mergers. We first generated two identical 
structures with Eddington's formula with /3 = 0. The two structures were collided, and the 
remnant of this collision was duplicated and collided again. This procedure was repeated 
until seven collisions were done. In the collisions the impact parameters were 0, and the 
collision axis changed from simulation to simulation^. The initial distances between the 
collided structures were 10.0, 12.5, 15.0, 25.0, 35.0, 40.0 and 50.0, respectively, and the 
structures had initial relative velocities of 0. 

There were 10^ dark matter particles in the first simulation and 1.28 • 10^ in the last 
simulation. We used units where the total mass and rg of the initial structure were defined 
to be 1, and a cutoff was made at a radius of 10. We ran the last simulation for 200 time 
units, so the structures had time to first collide (it happened after ~ 75 time units) and then 
form a new equilibrated halo. 

The simulations were done for initial structures following three different density profiles: 
a Hernquist profile, p{r) oc 1/(1 + r)^ and Eq. (2.3). 

Additional simulations were run with the same three density profiles, but instead of 
rotating the collision axis in the last simulation we used the same collision axis in simulation 
6 and 7^. 

2.6 Simulation V — Unstable Osipkov-Merritt models 

It is established that Osipkov-Merritt models following Hernquist profiles are unstable for 
low values of ran [37], due to the onset of a radial orbit instability [38], which creates a bar 
structure in the center of the halo. In this simulation we will study how such an instability 
affects the velocity anisotropy of a system. The perturbations in this simulation are therefore 
non-spherical because of the formation of a bar, despite the fact that initial conditions are 
spherically symmetric. 

We ran simulations with two density profiles; a Hernquist profile and Eq. (2.3). In both 
cases we defined rg = 1 and ran = 0.2rs. The total mass was 1. The two simulations were 
run for 200 time units and 10^ collisionless particles were used to represent each halo. 

2.7 Simulation VI — A halo and its subhaloes 

To see how cosmologically realistic haloes are affected by the dynamical friction from the 
subhaloes inside them, we generated a halo and subhaloes similar to what is found in the 
Via Lactea II simulation, which contains a main halo with a mass of 1.94 x IO^^Mq and a 
tidal radius of 462 kpc [39]. 

We used 10^ dark matter particles to represent the halo and the subhaloes, but only 
subhaloes heavier than 1.94 x 10® (the mass of 100 particles) were included. In total 7% 

'^The collision axes were the following in the seven simulations: x, y, z, x, y, z, x. The structures were not 
rotated between the simulations. 

^So the collisions occurred along the following axes in the mergers: x, y, z, x, y, z, z. 
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of the mass was contained in subhaloes. All haloes followed Hernquist density profiles with 
the same radius of maximum circular velocity, tidal radius, tidal mass, position and velocity 
as in the public catalog from the Via Lactea II simulation. The velocities of the main halo 
particles and the subhalo particles were sampled from velocity distributions calculated using 
Eddington's formula with /3 = 0. We ran the simulation for 100 time units, which is the same 
as 100 dynamical times for a particle at the scale radius. 

3 Analysis methods 

3.1 Analysing particles in cones 

In our analysis we will study particles in cones pointing in different directions. Each cone has 
an apex angles of 45°. For each particle and axis we calculate 9 = arccos(n • r/|r|), where n 
is a unit vector pointing in the direction of a cone, and r is a particle's position vector (seen 
from the center of the halo). Particles are selected to be contained in a cone, \i 6 < 22.5°. 
Analysing structures in this way was also done in [26]. 

We created 48 cones centered on the most bound particle of each halo. To define 
the pointing direction of each cone we generated 48 points on a sphere. The points were 
distributed in the same way as a sphere is divided into pixels in typical studies of the cosmic 
microwave background [40]. Note, that some cones will be overlapping each other in space, 
because of our choice of apex angle and number of cones. 

3.2 Determination of halo shapes 

Halo shapes and orientations are determined from the eigenvectors and eigenvalues of the 
moment of inertia tensor. First the shape in the central part is fitted with an ellipsoid, and 
next the shape is determined as function of radius with the procedure described in detail in 
[41]. In the remaining part of this article we will refer to the shape found with this method 
as the shape of a halo. In the following section we will also discuss contours of constant p 
and (3 in haloes, and in these cases we will explicitly mention, what kind of contour we are 
working with, so the different types of contours can be clearly distinguished. 

Note that we in this article let r denote the physical radius, and not the elliptical radius 
as it is done in some other studies. 

4 The behaviour of velocity anisotropy and halo shape 

To analyse the behaviour of shape and the direction-dependence of the velocity anisotropy for 
our haloes, we will present plots of /3(r) and p{r) through cones for a representative selection 
of haloes from our simulations. 

4.1 Category 1: Haloes with well-behaved /3-profiles 

Figure 1 shows halo H from simulation II, where structures were perturbed by varying G 
in time. Each of the lines represent a property (/3 or p) calculated through a cone for the 
structure. Also shown (in the right panel) are the axis ratios, c/a and 6/a, calculated from 
the inertia tensor. The plotted radii are in units of r_2 = r(7 = —2). 

The right panel shows the prolate shape of the halo. From the panel that shows p in 
different cones, it can also be seen that the halo is elongated since there is an offset between 
the densities in different cones. From the plot of (3 in different cones, it is seen that /3 is 
systematically larger along the major axis than along the minor axis at any given radius. 
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logr/r_2 log'r/r_2 log'r/r_2 



Figure 1. Velocity anisotropy profiles {left panel) and densities {central panel) calculated in 48 cones 
for halo H (a Hernquist profile with /^initial = 1— r^/(l + r^)) from the simulation, where G was 
changed instantaneously. Each line corresponds to a cone, and the coloured lines corresponds to 
cones which coincides with one of the principal axes. The right panel shows the axis ratios as function 
of radius. The cyan horizontal line segment in the left {central) panel shows the width of a contour 
with /3 = 0.5 (logp= -4.17). 

The contours of constant /3 are more elongated than the contours of constant p. The 
minor to major axis ratio of each contour can be approximated by the length of the cyan 
line segment in figure 1, and it is found that the /3-contour has a minor to major axis ratio 
of 0.1, and that the p-contour has a minor to major axis ratio of 0.6. 

Halo A from simulation II, where a sudden change of particle energies was made, is 
shown in figure 2. This halo is also close to being prolate, and we again see that the largest 
velocity anisotropy is in the direction of the major axis^. The minor to major axis ratio of 
the /3-contour is 0.35 and the axis ratio from the density shape is 0.49, so we again see that 
the /3-contour is more elongated than the p-contour. 

4.2 Category 2: Haloes with complex /3-profiles 

For halo H from the collapse simulation (figure 3) a slightly more complicated behaviour of 
/3 is seen. In the inner parts /3 is smallest along the minor axis of the density ellipsoid, and 
in the outer parts it is largest along this axis. Such a crossover is absent in the plot of p 
through different cones. Because of the crossover, the /3-contours of the halo do not have an 
ellipsoidal shape. 

Figure 4 and 5 show more examples of haloes with /3-profiles that are clearly not ellip- 
soidal. For the merger remnant /3 is roughly constant and positive along the major axis and 
monotonically increasing along the minor axis. This behaviour is consistent with our previ- 
ous study of major mergers [26]. The instability simulation gives a complicated behaviour of 
both /3, p and the axis ratios. Is is e.g. seen that p has significant wiggles and bumps along 
the minor axis. 

4.3 Category 3: The spherically symmetric haloes 

Finally we will also note that perturbed haloes also can follow simple spherically symmetric 
distribution functions. Figure 6 shows halo A from simulation I, and it has a spherical shape, 
and /3-profiles behaving similarly from cone to cone. 

^ When 48 cones with apex angles of 45° are distributed on a sphere, it can happen that the minor axis (or 
one of the other principal axes) is contained in several cones, since some of the cones are overlapping. This is 
why several cones are marked as being along the minor axis in Figure 2. 
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log r /r_2 log i'It-2 log t li'-2 



Figure 2. This plot shows properties of halo A (a Hernquist profile with Anitiai = 0) from simulation 
II. As in figure 1, a larger /? is found along the major axis compared to the minor axis at all radii. 




logr/r_2 logr/r_2 logr/r_2 



Figure 3. Halo I (a Hernquist profile) from the cold collapse simulation. In the inner parts /3 is 
largest along the major axis, whereas it is larger along the minor axis in the outer parts. 




Figure 4. Halo K (a Hernquist profile with initial = 0) from the merger simulation. Along the major 
axis /3 is roughly constant in the inner parts, and along the minor axis it is monotonically increasing. 



4.4 A summary of the categories 

To summarise the behaviour of velocity anisotropy and shape in our haloes, we have grouped 
them into three categories: 

1. One category is the triaxial haloes with /9 differing from cone to cone. The contours 
of constant (3 are typically more elongated than the constant-/? contours. Examples of 
such haloes are shown in figure 1 and 2. 
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logr/r_2 logJ'/r_2 logJ'/r-o 



Figure 5. Halo R(pcxl/(l + r)^) from the instability simulation. The density profile of this halo has 
several bumps and wiggles along the minor axis. /3 also exhibits different behaviour along different 
axes. 




logr/r_2 logr/r_2 logr/r_2 



Figure 6. Halo A (a Hernquist profile with /^initial = 0) from the G-variation simulation. This halo 
has a spherical shape and a spherically symmetric velocity anisotropy. 

2. We also see examples of triaxial haloes, where the constant-/? contours differ signifi- 
cantly from triaxial elhpsoids (see figure 3-5). We do for example see haloes, where /3 
is largest along one axis in the inner parts, and smallest along this same axis in the 
outer parts. 

3. In our simulations we also see almost spherical structures, where (3 is behaving similarly 
along each axis (see figure 6). 

Note, that the above list is not supposed to be a comprehensive list of all possible 
configurations of (3- and p-profiles in dark matter haloes. The purpose of the list is to 
highlight features, which have not been extensively studied in the literature. 

5 Relating the results to cosmological haloes 
5.1 Comparison with a cosmological halo 

A study which also examines halo shapes and the direction-dependence of /3-profiles is 
the analysis of the Via Lactea II halo, where the probability distribution of the velocity 
anisotropy, P{(3), is calculated in four different radial bins, and along the major, intermedi- 
ate and minor axes (see figure 7 from [27]). The Via Lactea II halo has a behaviour, which in 
several ways is consistent with our merger simulation in figure 4. In both studies (3 is largest 
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along the major axis (/3 ~ 0.3 — 0.4) and smallest along the minor axis (—0.6 ^ /? ^ 0.0) in 
the inner parts. The values of /? along the intermediate axis are somewhere in between. The 
halo from the Via Lactea II simulation is close to being prolate (c/a = 0.52 and h/a = 0.62), 
which is also similar to the behaviour of our merger remnant. Due to the similarities with our 
merger remnants, we conclude that the Via Lactea II halo is in category 2 of the categories 
from section 4.4. 

Merging is not the only process that can produce a halo, which is similar to the Via 
Lactea II halo. An example is the simulation of the unstable halo (figure 5), which has a 
/3-profile that behaves similarly to the mergers in the inner parts. The unstable halo is, 
however, significantly more elongated than the Via Lactea II halo. 

We encourage future studies of /3-profiles of cosmological haloes to use cones, since it 
reveals information about halo dynamics, which is hidden when properties are calculated in 
spherically averaged bins. 

5.2 Spatial anisotropy of galaxy kinematics in galaxy clusters 

In a recent study of galaxy clusters the line of sight velocity dispersion is calculated along 
the major and minor axis for a stacked sample of 1743 galaxies from the SDSS survey [42]. 
It is found that a significantly larger velocity dispersion is found for galaxies with positions 
along the major axis. The difference of the velocity dispersions along the two axes is 38 it 13 
km/s. A possible explanation for this difference is the direction-dependence of the velocity 
anisotropy profile, which we have studied in this article. A different explanation could be 
that the cluster's non-spherical mass-distribution can produce different line-of-sight velocity 
dispersions along the minor and the major axis, even though (5 has the same behaviour along 
these axes. 

6 Spherically averaged profiles 

In this section we will report the spherically averaged profiles of /3 and 7 together with the 
halo shapes. 

6.1 Velocity anisotropy profiles 

Different trends are seen in the spherically averaged velocity anisotropy profiles of our haloes, 
see figure 7. The setups in simulation I, II, III and V produce haloes with /3 ~ in the inner 
parts and /3 ~ 0.8 in the outer parts. This is consistent with what is e.g. reported in cold 
collapse simulations [36, 43]. 

The spherically averaged velocity anisotropy profiles of the merger remnants (simulation 
IV) are typically increasing from /3 = in the center of the halo out to the radius with (3 = 0.3, 
where a maximum appears. This behaviour is similar to many cosmological haloes [11], and 
in agreement with what is reported in other studies of merger remnants [19, 26, 44]. 

In the substructure simulation (simulation VI) only a tiny evolution of f3{r) is seen, 
which means that the dynamical friction from subhaloes only have a minor effect on velocity 
anisotropy profiles. This is consistent with the finding that the density profiles of cosmological 
haloes are not perturbed by the substructure [45]. 

6.2 7-profiles 

The six different physical processes in our simulations produce different 7(r)-profiles, see 
figure 8. In simulation I and II a core where 7 ~ appears in the inner regions (with 
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Figure 7. The spherically averaged velocity anisotropy profiles for the simulated haloes. 



logr/r„2 ^ —0.5). Such a core is absent in the remaining simulations. In the outer regions 
7(r) seems to be around 7 ~ —3 in simulation I, II and III. 

6.3 Halo shapes 

A huge diversity is seen in the halo shapes in figure 9. Several of the haloes from simulation 
I (with G-variations) are spherical and others are close to being prolate. In simulation II, 
III, IV and V the shapes of the haloes are close to being prolate. The halo in simulation VI 
(the simulation with subhaloes) remains spherical. 
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Figure 8. 7(r) for all the haloes. 
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Figure 9. The axis ratios, c/a {solid lines) and b/a {dashed lines), for the haloes, b/a+1 is plotted 
instead of b/a for presentation reasons. 
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Figure 10. Comparison between the attractor from HJS (i.e. simulation II) in grey^ and the other 
simulations. Wc have defined n = dloga"^,^^/ dlogr and 7 = dlogp/dlogr. In Table 2 the different 
colours and symbols can be found. Left panel: The structures exposed to the perturbations, where 
G was changed instantaneously. Central panel: The difference between the red lines and the blue 
lines is the orientation of the merger axis in the last simulation (sec section 2.5 for details). Two cold 
collapse simulations are also shown. Right panel: The two haloes from the instability simulation (V) 
and the subhalo simulation (VI). 

7 Comparison with the attractor 

An aim of the previous sections has been to examine the many different behaviours of halo 
shapes and /3-profiles of dark matter haloes. In the current section we will instead look for 
similarities of the different haloes, and compare them with the attractor from HJS. 

7.1 Motivation for the attractor 

The three quantities, /3(r), 7(r) = dlog p/(ilog r and K{r) = dloga"^^/ d\ogr, are clearly of 
importance for spherical and static distributions of collisionless structures, since the Jeans 
equation [32] can be written as 



This motivates an analysis of structures behaviour in the (/3, 7, K)-space. In the HJS-paper 
such an analysis was made for the structures from simulation II (the energy exchange simu- 
lation), and in this 3-dimensionaI space all haloes followed a 1-dimensional relation, i.e. an 
attractor. In this section we will compare this attractor prediction with our structures. 

7.2 Comparing haloes v^ith the attractor 

First we will analyse structures in the (/3, 7 + K)-projection, see figure 10. The points show 
properties calculated in spherical bins distributed logarithmically in radius. For simulation 
I, where G was changed multiple times, the different haloes end up on a one-dimensional 
curve, like in HJS. The scatter at large /3 is smaller than in the original HJS-simuIations, 
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and there is an offset towards smaller (7 + K)-values in simulation I compared to the HJS- 
simulations (simulation II) . The same trend is seen in the study of the attractor by [46] . The 
two haloes from the cold collapse simulations (simulation III), also end up on the attractor 
in this projection. 

Simulation IV involves multiple major mergers. The result of the three simulations 
with different density profiles are shown as red points in figure 10. The blue points show 
the additional simulations with identical collision axes in the last two simulations. The 
spherically averaged profiles follow the attractor in the inner parts, even though /3 behaves 
different along each axis (see figure 4). When inspected the attractor is not obeyed when /3, 
7 and K are calculated in cones. 

The unstable model in Simulation V ends up close to the attractor. The halo in simula- 
tion VI is not perturbed enough by the presence of the substructures to be dragged towards 
the attractor. 

7.2.1 Different projections 

In figure 11-13 the structures are plotted in the three principal projections, i.e. (7, /3), {k,/3) 
and (7, k). In all projections there is a very good agreement between simulation I and II, so 
from now we will define the attractor by the result of these two simulations. 

In the (7, /3)-projection (figure 11) the cold collapse process produces haloes, which are 
consistent with the attractor prediction. The spherically averaged properties of the merger 
remnants follow the attractor in the inner parts. We also see that the haloes from simulation 
I-IV follow the linear 7-/3 relation from [24] in the inner parts, where 7 > —2.2. The unstable 
model is not on the attractor. 

In the (k, /3)-projection (figure 12) the collapse simulations are again close to the attrac- 
tor prediction, but some significant deviations are present in the inner parts with /3 < 0.2. 
The merger simulations behave differently and the unstable models only follow the attractor 
(in this projection) in the outer parts with /3 > 0.35. 

7.2.2 Comparison vi^ith cosmological pseudo-phase-space densities 

Figure 13 shows the simulations projected onto the (7, K)-plane, together with the relation, 
K = |(7 -|- q), which follows from the connection, 

p/aUr) oc r-, (7.2) 

with a = 1.91 [11], which comes from fitting cosmological haloes over their entire radial 
ranges. 

We see that the attractor prediction is inconsistent with the relation from Eq. (7.2), 
which is obeyed by cosmological haloes. It can also bee seen that the inner parts of the 
merger remnants and the haloes from the collapse simulations are closer to relation (7.2) 
than the attractor prediction. 

7.2.3 Energy exchange in the relaxation processes 

The purpose of this subsection is to study the change of particle energies in the relaxation 
processes. For our structures we identified the particles, that ended up in a thin spherical 
shell centered at r_2 = ^(7 = —2) in the last snapshot, and monitored how their energies 
changed through the simulation. In figure 14 we have plotted the dimensionless quantity 
*dyn(''~2) X {\{dEi / dt) / Ei\) [Ei \s the energy of the z'th particle, and (• • • ) denotes the mean 
of the particles that ended up in the shell at r_2) as function of t/idyn('^-2)- 
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Figure 11. A projection in the (7, /3)-planc. Also plotted is the linear relation {thick green line) from 
[24]. 
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Figure 12. The (k, /3)-projection. 

For the cold collapse experiment (the structure following a Hernquist profile is shown), 
the collapse begins immediately after the start of the simulation, and a redistribution of the 
energy occurs. For the merger experiment (the Hernquist profile is shown) the exchange of 
energy grows in the beginning, where the two haloes approach each other, and the largest 
peak occurs approximately when the two haloes collide for the first time {t ~ 40tdyn)i ^-nd 
a smaller peak is present when they collide for the second time {t ~ 55tdyn)- At later times 
the energy exchange is more than an order of magnitude lower than the peak values for both 
experiments. 

The energy exchange in the simulation with an unstable model (the structure following 
Eq. (2.3) is shown) is weaker than in the collapse simulation and the mergers. The instability 
causes the energy-exchange around the peak at t = 22tdyn- When inspecting the /3-profile in 
each individual snapshot it is seen, that it changes shape from t/tdyn = 16 until i/^dyn = 24. 
In the later snapshots no significant evolution of this halo is seen. 
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Figure 13. The (7, K)-projection. 




Figure 14. The fractional change in energy per dynamical time for particles that end up in a spherical 
shell around r_2 in the final snapshots. In the cold collapse simulation and the merger simulations 
large changes in the energy of the particles appear. Much lower variations in the energy are present 
in the simulation with the substructure. For reference an equilibrium structure with no substructure 
is also shown. 

In the subhalo simulation a small energy exchange is present throughout the simulation. 
To see the amount of energy exchange generated by the substructure, we ran a similar 
simulation with a halo in equilibrium; i.e. a halo where the subhaloes were removed, and 
ordinary particles were distributed throughout the halo to compensate for the removed mass. 
By comparing the two simulations we find that only the small bump at t < 5tdyn; where some 
of the particles in the substructure are stripped, might be caused by the subhaloes. 
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7.3 An overview: which structures are on the attractor? 

The overall pattern in our simulations is that the haloes that end exactly on an attractor 
are those from simulation I and II, which involved artificial processes that changed the grav- 
itational constant or the kinetic energy of the particles. The collapse simulations ended up 
close to the attractor, but not exactly on it, and the merger simulations had clear deviations 
in the outer parts. 

We also found that haloes following the attractor deviate from the universality, p/cr^^^ oc 
r~", which is found in cosmological simulations. Our conclusion is that the dynamics and 
structure of cosmological haloes are different from the attractor. Further indications, that 
cosmological haloes are not on the attractor, comes from the finding that our cold collapse 
and merger simulations, which are both cosmologically realistic processes, produce haloes 
that are not exactly on the attractor. The strong direction-dependence of (3 found in the Via 
Lactea II halo [27] also supports that cosmological haloes are not on the attractor. 

For both the merger simulations and the collapse simulations we do, however, see that 
the /3- and 7-profiles (spherically averaged) in the inner parts of the structures obey the at- 
tractor in the (7, /3)-projection. Since a very good agreement only is found in this projection, 
cosmological haloes are better described by the 7-/3 relation from [23] than by the attractor. 

It is important to note that this 7-/3 relation only applies to the spherically averaged 
values of (3- and 7-profiles. The haloes still have freedom to have different /3-profiles in 
diff'erent directions as long as the spherically averaged profiles obey this relation. An example 
of a halo in which the /3-7 relation is not obeyed through cones in different directions, but 
only in spherically averaged bins, has been presented in the study of merger remnants in [26]. 

7.4 The relevance of the attractor 

We have shown that typical cosmological haloes have departures from the attractor predic- 
tion. Even though this is the case, we still believe that there must be a physical origin of 
the similar behaviour of the haloes from simulation I and II, when they are analysed in the 
(/3, 7, K)-space. Understanding this behaviour could potentially lead to a better understand- 
ing of violent relaxation and mixing processes in collisionless systems. 

8 Summary 

In this article we have examined velocity anisotropy profiles and halo shapes in a range of 
pure dark matter simulations. We have found that haloes with elongated shapes can have 
several kinds of direction-dependent velocity anisotropics. In some cases the /3-profile is 
aligned with the density profile, and the largest velocity anisotropy is found along the major 
axis. In other cases the /3-profiles are more complicated, and it is not possible to define an 
axis along which the velocity anisotropy is largest at all radii. Such a behaviour is e.g. seen 
in the remnant of a major merger, which has a shape and a direction-dependence of /3, that is 
very similar to the Via Lactea II halo. We suggest that future studies of cosmological haloes 
should calculate the velocity anisotropy profiles in cones, since it reveals properties of the 
true velocity distributions, which are hidden from the spherically averaged profile. 

We have also compared our haloes with an attractor [25], and we conclude that cos- 
mological haloes have departures from the attractor prediction. We do, however, find that 
the spherically averaged (3- and the 7-profile (7 is the logarithmic derivative of the density 
profile) obey an approximately linear relation (from [23]) in the inner parts, even though 
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this relation is not necessarily obeyed for particles in individual cones pointing in different 
directions. 
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